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Abstract. Modulational instability in a photorefractive medium is studied in the 
presence of two wave mixing. We then propose and derive a model for forward four 
wave mixing in the photorefractive medium and investigate the modulational instability 
induced by four wave mixing effects. By using the standard linear stability analysis 
the instability gain is obtained. In both the cases, the geometry is such that the effect 
of self-phase-modulation self-focusing is suppressed and only the holographic focusing 
nonlinearity is acting. 
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1. Introduction 

A soliton is a localized wave that propagates without any change through a nonlinear 
medium. Such a localized wave forms when the dispersion or diffraction associated 
with the finite size of the wave is balanced by the nonlinear change of the properties of 
the medium induced by the wave itself pQ. Solitons are ubiquitous in nature and have 
been identified in various physical systems like fluids, plasmas, solids, matter waves and 
classical field theory. Such self-trapped beams are the building blocks in future ultra- 
fast all-optical devices. They can be used to create reconfigurable optical circuits that 
guide other light signals. Circuits with complex functionality and all-optical switching 
or processing can then be achieved through the evolution and interaction of one or more 
solitons [2]. It is possible to compress circuitry into a compact space with many circuits 
sharing the same physical location. Furthermore, certain photosensitive materials offer 
the potential for erasing one light written device and replacing it by another. Hence, 
we have the building blocks for dense reconfigurable virtual circuitry. Due to the 
development of materials with stronger nonlinearities the optical power needed to create 
such virtual circuits has been reduced to the milliwatt and even microwatt level, bringing 
the concept nearer to practical implementation. In a recent experimental study [3], 
optical spatial screening solitons have been observed for the first time in strontium 
barium niobate (SBN). Soliton states of this sort are known to occur when the process 
of diffraction is exactly balanced by light-induced photorefractive (PR) waveguiding [I]. 

Modulational instability (MI) is a process in which tiny phase and amplitude 
perturbations that are always present in a wide input beam grow exponentially during 
propagation under the interplay between diffraction (in spatial domain) or dispersion 
(in temporal domain) and nonlinearity. Instabilities and chaos can occur in many 
types of nonlinear physical systems. Optical instabilities can be classified as temporal 
and spatial instabilities depending on whether the electromagnetic wave is modulated 
temporally or spatially after it passes through the medium. Temporal instability has 
been studied by various authors El [7] , and the first experimental observation of MI in 
a dielectric material was in 1986 [8]. The temporal MI occurs as an interplay between 
self-phase modulation and group velocity dispersion. In optical fibers, MI occurs in 
the negative dispersion region and is responsible for the breakup into solitons. In the 
spatial domain, diffraction plays the role of dispersion. When both diffraction and 
dispersion are present simultaneously, it results in spatio-temproral MI [9j. Recently, 
Wen et al. investigated MI in negative refractive index materials |10| . In a rather loose 
context, MI can be considered as a precursor of self-trapped beam formation. In the 
spatial domain, MI manifests itself as filamentation of a broad optical beam through 
the spontaneous growth of spatial-frequency sidebands. The MI is a destabilization 
mechanism for plane waves. It leads to derealization in momentum space and, in turn, 
to localization in position space and the formation of self-trapped structures. During 
MI, small amplitude and phase perturbations tend to grow exponentially as a result of 
the combined effects of nonlinearity and diffraction. This results in the disintegration of 
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a large-diameter optical beam during propagation. Such instabilities have been widely 
studied in photorefractive media [HJ [121 EE] ■ Assanto et al. [H] investigated transverse 
MI in undoped nematic liquid crystals, a highly nonlocal material system encompassing 
a reorient at ional nonlinear response. They observed the one- dimensional development of 
transverse patterns which eventually lead to beam breakup, filamentation, and spatial 
soliton formation, using both spatially coherent and partially incoherent excitations. 
The MI phenomena has been previously observed in various media like Kerr media 
[13] . electrical circuits [16J, plasmas [17] . parametric band gap systems [3J, quasi-phase- 
matching gratings [18], discrete dissipative systems [19] and in PR polymers [20]. The 
transverse instability of counterpropagating waves in PR media is studied by Saffman 
et al [21] . From the above investigations, it is clear that the study of MI in a medium 
is both of fundamental as well as of technological importance. 

In this work, we investigate the modulational instability induced in a 
photorefractive medium by two wave mixing and the forward four wave mixing occurring 
in the PR material. The photorefractive medium is a well suited medium for the 
practical implementation of virtual circuits as self-trapping of beams can be observed 
in this medium at very low laser powers. Four wave mixing in the phase conjugate 
geometry is widely used in photorefractive materials for various applications. It can be 
used to implement several different computing functions [22J, optical interconnects |23j, 
matrix addition [21], and optical correlator [25]. Nonlinear solutions for photorefractive 
vectorial two-beam coupling and for forward phase conjugation in photorefractive 
crystals have been found in [26]. Recently Jia et al. [27] experimentally demonstrated 
degenerate forward four wave mixing effects in a self-defocusing PR medium, in both one 
and two transverse dimensions. They observed the nonlinear evolution of new modes as 
a function of propagation distance, in both the near-field and far-field (Fourier space) 
regions. 

The motivation behind studying MI in such a geometry is due to the recent proposal 
by Cohen et al. [28] of a new kind of spatial solitons, known as the holographic (HL) 
solitons. They are formed when the broadening tendency of diffraction is balanced by 
phase modulation that is due to Bragg diffraction from the induced grating. Holographic 
solitons are solely supported by cross-phase modulation arising from the induced 
grating, not involving self-phase modulation at all. In 2006 [29], they showed that 
the nonlinearity in periodically poled photovoltaic photorefractives can be solely of the 
cross phase modulation type. The effects of self-phase modulation and asymmetric 
energy exchange, which exist in homogeneously poled photovoltaic photorefractives, 
can be considerably suppressed by the periodic poling. They demonstrated numerically 
that periodically poled photovoltaic photorefractives can support Thirring-type (solitons 
which exist only by virtue of cross phase modulation) (holographic) solitons. HL solitons 
in PR dissipative medium was studied by Liu [30J. Existence of HL solitons in a 
grating mediated waveguide was studied by Freedman et al. [31]. Salgueiro et al. 
[32] studied the composite spatial solitons supported by mutual beam focusing in a 
Kerr-like nonlinear medium in the absence of the self-action effects. They predicted the 
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existence of continuous families of single and two-hump composite solitons, and studied 
their stability and interaction. In 2007, two-dimensional holographic photovoltaic bright 
spatial solitons were observed in a Cu:K .25Nao.75Sri. 5 Bao.54Nb50i5 crystal in which two 
coherent laser beams, a signal beam, as well as a strong and uniform pump beam at 532 
nm are coupled to each other via two- wave mixing [33] . 

This paper is organized as follows. The basic propagation equation for the two 
wave mixing geometry is presented in Section [2] and the modulational instability in 
this geometry is studied. In section [3] the governing equations for the forward four wave 
mixing geometry is presented. The system is studied without using the undepelted pump 
beam approximation. The standard linear stability analysis for the coupled equations 
is carried out and the gain spectrum is obtained. Section 0] concludes the paper. 



2. Two wave mixing geometry 

In this section, we present the derivation of the model equation for the two wave mixing 
geometry [3U [35] and study the MI in it. Consider the interaction of two laser beams 
inside a photorefractive medium. Stationary interference pattern is formed, if the two 
beams are of the same frequency. 

Let the electric field of the two beams be written as 

Ej = Aj exp[i(ujt - kf.r)], for j = 1, 2. (1) 

Here, A\ andv4 2 are the amplitudes, ui is the angular frequency, and ki and k 2 are the 
wave vectors. 

The medium is assumed to be isotropic and both beams are polarized perpendicular 
to the plane of incidence. The total intensity of the beams 

/= \E\ 2 = \E 1 + E 2 \ 2 , (2) 

can be expressed as 

1= \A 1 \ 2 + \A 2 \ 2 + AiA£exp[zK.r] + A 2 A[exp[-iK.r], (3) 

where 

K = k 2 — ki. 

The magnitude of the vector K is 27r/A where A is the period of the fringe pattern. 

Equation [3] represents a spatial variation of optical energy in the photorefractive 
medium. Such an intensity pattern will generate and redistribute charge carriers. As a 
result, a space charge field is created in the medium. This field induces an index volume 
grating via the Pockels effect. In general, the index grating will have a spatial phase 
shift relative to the interference pattern. 

The index of refraction including the fundamental component of the intensity- 
induced gratings can be written as 

'n x A\A 2 



n = n + 



2/n 



■ exp[i0] exp[— zK.r] + cc 



(4) 
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where no is the index of refraction when no light is present, is real and rii is a real 
and positive number. The wave equation reduces to the Helmholtz equation for highly 
monochromatic waves like laser beams which can be viewed as a superposition of many 
monochromatic plane waves with almost identical wave vectors 

V 2 E + oo 2 n 2 E/c 2 = 0, (5) 

where E = E\ + E 2 . 
Let 

Ej = Aj(x, y, z) exp[iuit — ifyz], 

where A(x, y, z) is the complex amplitude that depends on position, p 1 and /3 2 are 
the z component of the wave vectors ki and k 2 inside the medium respectively and z is 
measured along the central direction of propagation. We solve for the steady state so 
that Aj is taken to be time independent, so that equation (jSJ) becomes 

d 2 d 2 d 2 \ 

+ L + l, + - 2 n 2 /c 2 )E = 0. (6) 



, dx 2 dy 2 dz 2 

Substituting equation PJ in equation fl§]) and solving by neglecting the second order term 
n 2 and using the fact that for highly directional monochromatic waves with a » A, 
d 2 A/dz 2 can be neglected, we get 

= Vl^i + ^le^A^A, (7) 

OZ C iq 

and 

= V\A 2 + ^le^Al^A,. (8) 

OZ C 2 I 

For the case when the two laser beams enter the medium from the same side at 

Pi — 02 — kcosO = —n cos6. 

A 

where, 28 is the angle between the beams inside the medium. Simplifying we obtain 

oz 2k cos 2I 
^ = ^ L -nVlA 2 + ^-A 2 \Ai\ 2 , (9b) 

oz 2k COS V 2in 



where 



2 ™i -i<f> 



. (10) 
A cos# 

is the complex coupling coefficient. A similar model was used to study the existence 
of Holographic solitons [28]. A prerequisite for obtaining Holographic focusing is that 
the induced grating be in phase with the intensity grating. If the induced grating is 
shifted by tt, then the grating leads to holographic defocusing. If the grating is ±7r/2 
phase shifted with respect to the intensity grating, then the interaction will yield an 
asymmetric energy exchange between the two beams because of the two beam coupling 
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property of the photorefractive material. Therefore, = for the existence of bright 
solitons. This gives the coupling constant as 

r =^V <»> 

A cos# 

The next step is to study the propagation of a broad optical beam through a PR 
medium. We study MI of a one dimensional broad optical beam. Hence the y dependent 
term in the transverse Laplacian in equation (19 a\) can be neglected. For a broad optical 
beam, the diffraction term in equation ( |9al ) can be set to zero giving the steady state 
solutions as 

r r? i 

Aj = VPexp - i—z , (12) 

L Zl J 

for j = 1,2. 

To study MI, we consider small perturbations of the steady state solutions : 

r TP i 

Aj = (vP + a,,) exp — i— -z . (13) 
L 2Iq . 

Substituting in (l9~a|) and neglecting the quadratic and higher order terms in aj, the 
perturbations a x and a 2 are found to satisfy the following linearized set of two coupled 
equations : 

.da x 1 d 2 a x TP 



"dz 2kcos6dxi + 2T {a2 + a ^ (14) 
.da 2 1 d 2 a 2 TP 

^ = 2^^ + 27o (ai + fll) - (15) 
It is important to note that the evolution of the perturbations depend solely on 
the cross phase modulation. To solve equations (TH|) and f|T5l) . we assume that the 
perturbations be composed of two side bands : 

aj(x, z) = Uj(z) exp[iKx) + Vj(z) exp[— inx]. (16) 

The substitution of equation (TIB"]) in (TTJJ and ffl5l) results in a set of four 
homogeneous equations in U±, U2, V% and V 2 as 

oz 2/ccost^ 21 Q 
8V* k 2 TP 

• ^ ./.-fe+yr), (19, 



dz 2k cos6 2I 
dV 2 * . k 2 „ .TP 



oz 2k cos 9 2I 

This set has a nontrivial solution only if the determinant of the coefficient matrix 

vanishes. The eigenvalues of the system are obtained as 

1 {-I k 4 -2ATPk 2 \Vz , , 

A 1± = ±— — °- , (21) 



2A\ I 
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where A = kcosd. The plane wave solution is stable if perturbations at any wave 
number k do not grow with propagation. MI gain will exist only when i?e[A] > 0. This 
condition is satisfied only by the second set of eigenvalues. The gain associated with 
the system is given by 

G = |^(A±)|. (23) 

A typical plot of the gain spectrum is given in figure ([1]). We consider the case of 
BaTiC>3 crystal with a space charge field of 10 4 V/m, electro-optic coefficient r^ = 
1640 x 10 -12 m/V and refractive index n = 2.4 which gives a coupling constant T = 20 
cm -1 . Figure [2] gives the variation of gain with the angle 9. The gain increases with 
decrease in angle. 

3. Forward four wave mixing geometry 

In the two wave mixing case, two coherent beams interfere inside a photorefractive 
medium and produce a volume index grating. In the case of optical phase conjugation, 
using four-wave mixing, a third beam is incident at the Bragg angle from the opposite 
side and a fourth beam is generated. We now consider the situation in which the third 
beam is incident from the front at the Bragg angle and a diffracted beam is generated. 
That is we consider the interaction of four beams in a PR medium in the forward 
geometry (see figure (jHJ)). We assume that all the beams have the same frequency to. 
We propose a model for the observation of MI in a PR medium induced by four wave 
mixing. The method proceeds by first writing the four coupled equations for the present 
geometry. Of the four beams, let beams A 2 and A3 be the pump beam, A± be the signal 
beam and A4 be the generated beam. Beam A\ is coherent with beam A 2 and beam 




Figure 2. Variation of gain with respect to the angle between the two beams. The 
lowermost curve is for 9 = n/16, the dotted curve is for n/8 and the topmost curve is 
for 7r/4. 




Figure 3. Forward four wave mixing in photorefractive media in the transmission 
geometry. 
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A 3 is coherent with beam A±. Thus the index grating consists of two contributions : 
A\A 2 and y4gy4 4 . The index of refraction including the fundamental component of the 
intensity-induced gratings can thus be written as 

"ni(AJA 2 + A*A±) 



n 



n + 



2/ 



exp[i0] exp[— iK.r] + cc 



(24) 



where k 2 — ki = k4 — k3 = K. 

Following the above procedure for the two wave mixing case, we get the following 
four coupled equations for the present geometry: 

.dA x 



dz 
.dA 2 

dz 
.dA* 

dz 
.dA± 



_L_ vlAl + I_ {AlAt + A3AdA2 , 

T^—K^lA, + ^{A lA \ + A 3A \) Ai , 
2k cos v 21 

^^IM + ^-{a\a 2 + a* 3 a,)a 3 , 

2k cos d 2io 



(25 a) 
(256) 
(25 c) 
(25 d) 



dz 

where I = h + I 2 + h + h- 

The model permits plane wave solutions of the form Aj(x,z) = \/Pexp[ 
The next step is to carry out a Linear Stability Analysis of the plane wave solutions. 
For this the plane wave solution is perturbed as 



•rp 1 



An 



P + aj(x, z)) exp 



(26) 



where a,j is a small complex perturbation. Inserting this into the coupled equation 
and linearizing around the solution yields the equations for the perturbations: 



.da\ 


r 


dz 




.da 2 


r 


dz 




.da 3 


r 


dz 


~2lo 


da± 


r 


dz 


~27o 



ai + 



a 2 + 



04 + 





1 


d 2 a x 


2k 


cos 


9 dx 2 




1 


d 2 a 2 


2k 


COS 


9 dx 2 




1 


d 2 a 3 


2k 


COS 


9 dx 2 




1 


<9 2 04 


2k 


COS 


9 dx 2 



+ 77p(2a 2 + a 2 + a 3 + a 4 ), 
2±o 



r 

2ir 



(2ai + a\ + ag + a 4 ), 



— (ai + a 2 
zi 



2a4 + a! 



■{)■ 



+ 7Tr( a i + a 2 + 2a 3 + ag). 

ZIq 



(27) 
(28) 
(29) 
(30) 



Now, we assume that the spatial perturbation a(x, z) is composed of two side band 
plane waves, i.e 

a,j(x, z) = Uj(z) exp[iKx) + Vj(z) exp[— ikx]. (31) 
Substituting, we get eight homogeneous equations as 

r 2 + V* + U 3 + V;), (32a) 
\ + V* + V; + U A ), (32b) 



dU x ., 


TP 


K 2 


)U X 


TP 


dz = l{ ' 


2I 


2k cos 9 




TP 


H 2 


)U 2 


.TP 
~ % 2h 


dz ~ l{ ' 


2/ 


2k cos 9 
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dz 



fiTL 




r p 


t-2 

IX 


ft "7 






9k rat f) 


f)TJ t 
(JU4 


—i(- 


r p 
1 1 


K 2 
IX 


it v 








u 


—i( 


V P 


K 2 


f) y 


2 In 


2 A* 00s 


dV* 


—i( 


TP 


/t 2 


dz 


2/ 


2fc cos 


dV 3 * 


— z( 


TP 


/t 2 


dz 


2/ 


2r€ COS 


dV 4 * __ 


1 1 


rp 


/t 2 



rp 
rp 



+ 2U 4 + V?) 



. l c/4-&(v 1 * + [/ 2 + 2c/ 3 + y 3 *) : 

)V* + z^(2V 2 * + U 2 + V 3 * + U 4 ) 

)v 2 * + i^(2v 1 * + u l + u 3 + v;) 



TP 
TP 



u 2 + 2v; + u 4 ), 



2J n 2A; cos 



)^4* - i^riUi + V* + 2V; + U 3 

Zln 



(32 c) 
(32 d) 
(32 e) 
(32/) 
(32<?) 
(32/i) 



The obtained eight coupled equations can be written in the compact matrix form 



*1T 



as, d z X = MX, where M is an 8x8 matrix with X = [U x U 2 U 3 U 4 V{ V 2 * V 3 * V 4 

This system has a nontrivial solution only if the determinant of the matrix vanishes. 
The real part of the eigenvalues of the stability matrix of equation ([3]) gives the gain 
associated with the system [36J. Out of the eight roots of the system, only the root with 
maximum positive value contributes to the MI gain. Here only the eigenvalue given by 



V=gg+4grT^gp 

A± = ± ■= (33a) 

2 V Ink cos 6 

contributes to the MI of the system. 

We note that the gain vanishes for all values of k greater than n max = AkTPj Jo- 
Defining 7 = TP I In, we can rewrite the gain as 



" " " ^'^-{—) )- (34) 



7 

The variation of the gain coefficient in the forward four wave mixing process 
with respect to the spatial perturbation is plotted in figure PJ. Such instabilities are 
useful for pattern formation. A transverse modulation instability of a single beam or 
counterpropagating beams is a general mechanism that leads to pattern formation in 
nonlinear optics [371 EH]- We expect that similar results will be obtained using the 
present geometry. 



4. Conclusion 



We first studied MI occurring in a rR medium in the two wave mixing geometry and 
further modeled the forward four wave mixing occurring in the rR medium and studied 
MI in this geometry. In both the cases, the geometry is such that only the holographic 
focusing nonlinearity is acting. MI does not rely on self-phase-modulation self-focusing 
but results only by virtue of the competition between induced periodic modulation of 
the refractive index and diffraction of the beam. The MI gain spectrum is obtained for 
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Figure 4. A typical plot showing the gain spectrum of the system in the forward four 
wave mixing process with respect to the spatial perturbation n. 

both two wave mixing and forward four wave mixing geometry. Such instabilities will 
be useful for pattern formation. Photorefractive materials are attractive for the studies 
of pattern formation as their slow time constant gives the possibility of observing the 
spatiotemporal dynamics of the system in real time. This also reduces the demands on 
experimental equipment where speed is often a crucial parameter. Another advantage of 
photorefractive pattern formation is that patterns can be observed using optical powers 
of tens of mW. In contrast, pattern formation using other nonlinearities require optical 
powers in the order of 1W. 
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